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Abstract: We point out that there is a parameter region in supersymmetry with heavy 
scalars and higgsinos, in which the heavier of bino and wino becomes long-lived as a conse¬ 
quence of the heavy higgsinos. In this region these electroweak gaugino sectors are secluded 
from each other with very small mixings that are inversely proportional to the higgsino 
mass, jji. We revisit the bino and bino decays and provide simple formulae for the partial 
decay rates and the lifetimes in the limit of heavy higgsinos. The scale of fi required for 
the long-lived electroweak gauginos highly depends on the mass hierarchy between bino 
and wino. The neutral wino can be long-lived (cr > 1 cm) even with |/x| ^10 TeV if 
~ 20 GeV. We discuss the collider signatures of the long-lived binos and winos 
in this scenario. 
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1 Introduction 

Supersymmetry (SUSY) remains as a promising new physics candidate after LHC Run-I. 
The discovery of the Higgs boson with nih — 125 GeV [1-4] and the negative results of 
the new physics searches at the LHC make the SUSY scenario with heavy scalars more 
attractive. In such a scenario, scalars (except for the Standard Model (SM) like Higgs 
boson, h) are heavier than the TeV scale and the observed Higgs mass can be easily realised 
by the large corrections from the heavy scalars [5-7]. The gauge coupling unification can 
be achieved with the light gauginos [8] and such light gauginos can be within the reach 
of the collider experiments. If the wino is the lightest SUSY particle (LSP), the thermal 
relic abundance bounds the wino mass from above by 2.7 — 3 TeV [9-11]. Phenomenology 
of light neutralinos with the thermal relic density smaller than the observed dark matter 
abundance was extensively studied in the literature [12-24]. The splitting between the 
scalar and gaugino masses is also motivated from model building perspective. If the SUSY 
breaking field with non-vanishing T-term is charged under some symmetry, only scalars 
acquire the soft masses at tree-level. The gaugino masses can be generated by other 
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mechanisms such as the anomaly mediation, but they are parametrically smaller than the 
scalar masses. 

It has been pointed out that gluinos can be long-lived in this scenario if squarks are 
heavier than 0(10^) TeV [ 8 ]. Such long-lived gluinos, if produced, may leave a distinctive 
i?-hadron signature or a signature of displaced vertices [25, 26]. With decoupled squarks, 
the SU (3) gaugino sector is secluded from the other gaugino sectors and the gluino becomes 
stable even if there exist lighter gauginos. The situation is somewhat different for the 
SU (2 )L and U{1)y gauginos because of the electroweak (EW) symmetry breaking (EWSB). 
These sectors are mixed after the EWSB and the heavier gaugino can decay into lighter 
ones without mediation of scalars. However, this mixing is proportional to the ratios of 
the vacuum expectation values of the Higgs field and the mass of the fermionic component 
of the Higgs field: Vu/^J■ and Vd/^^, and can become arbitrarily small if the magnitude of jjL 
increases. 

In this paper we point out that there is a region in SUSY parameter space with 
heavy scalars and higgsinos where the heavier of EW gauginos becomes long-lived as a 
consequence of the large ^ parameter. We revisit the decays of EW gauginos and derive 
simple formulae for partial decay rates and lifetimes, which are valid in the limit of the 
heavy higgsinos. We identify the region where the bino and winos become long-lived for 
both cases of two- and three-body decays. 

The paper is organised as follows. In the following section we review the chargino and 
neutralino sectors of the minimal supersymmetric Standard Model (MSSM) and study the 
mixings and interactions between bino and winos in the heavy higgsino limit. We then 
revisit the two- and three-body decays of bino and winos in sections 3 and 4, respectively. 
In these sections, we identify the scale of /U, at which the bino and winos become long- 
lived in collider experiments. In section 5 we briefly discuss how the large mass splitting 
between the gauginos and higgsinos can be theoretically achieved. The collider signatures 
of the long-lived bino and winos in our scenario is briefly discussed in section 6 . Finally, 
we conclude in section 7. In appendix A we list detailed formulae for auxiliary functions 
used in the paper and in appendix B we compare approximate expressions obtained here 
with an exact calculation using SDecay [27]. 


2 Interactions between winos and bino 

The chargino mass matrix is given by;^ 


( M 2 \f2mws^ 

■^C = n. 

yV2mwCi3 /r 

where C /3 (s^) represents cos [3 (sin /3) and /? is defined as the ratio of the vacuum expectation 
values of two Higgs fields, tan/3 = {Hu)/{Hd). In the large /r limit, this matrix can be 

^ We closely follow the calculation and conventions used in [28, 29]. 
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diagonalised as = U*AicV ^ where U and V take the following forms 


1 


\/2sjsnrnY 

[/ ~ I ^ ^ 

\/ 2 spmyv ^ 


V ~ 



The eigenvalues are then given by: 

^W± - ^2, 

with m^± <C 

The neutralino mass matrix is given by 


m-± ~ |//|, 


( 2 . 2 ) 


(2.3) 


Mjv = 


/Ml 0 

0 M 2 

-mzswC/3 ruzcwcp 


-mzswci3 mzswsp \ 
mzcwcp —mzcwSji 
0 —^ 


\ mzswSjS -ruzcwsp 0 

In the large [i limit, the diagonalising matrix Z can be written as 

( 


(2.4) 


/ 


Z ~ 


1 

'm%S2WS2p 

mzswS/S 

2/^(M2—Ml) 

h 


1 

rrizcwsp 

2fi{M2—Mi) 

h 

mzsw(cp—sp) 

razcw{cp—sp) 

1 

\/2/i 

V^/r 

V2 

mzswicp+sp) 

mzcw{cp+sp) 

1 

\/2/i 


V2 


1 

\/2 

1 

x/2 


(2.5) 


/ 


where sw is the sine of Weinberg angle, sin 6 * 14 /) and cw = cos^vv- The mass eigenvalues 
are found to be 

~ Ml, m^o ~ M 2 , m^o ~ m^o ~ |/i| , (2.6) 

where <C m^o,m^o. 

The wino and bino interact with each other and a Higgs or a gauge boson through the 
mixing matrices U, V and Z. The interaction with the light Higgs boson, h, is dictated by 
the X^-Xrh coupling [29] 


G L _ _ 

ijh — “ 


ijh 


2 sw 


{Zj 2 — Zji tan 6 \Y){Zi 3 sin a + cos a) + i j , 


(2.7) 


where a is the mixing angle of the neutral Higgs sector and reduces to sin a = cp, cos o. = 
—in the decoupling limit of the SUSY Higgses [30-32]. In the large [i limit these couplings 
can be written as 


G 


W°Bh 


T, - ( c / 3 - 2’13 - SpZii) - —{cpZ23 - SpZ2A) 

Isw L Cw 


'mzS 2 p 


( 2 . 8 ) 


This interaction is illustrated diagrammatically in the top-row graphs of figure 1. We 
note that the coupling originates from the bino-higgsino and wino-higgsino mixing and is 
suppressed by mz/fJ- for large 

^In order to stress wino nature of the lightest chargino, we denote it here as rather then usual xt- 
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Figure 1. The wino-bino interaction with the SM Higgs and gauge bosons. 


The interaction with the W boson is described by the coupling 

- —Zi4Vj2), G^w = — {ZaUji + —^ZaUj2) , 

sw v2 sw v2 


(2.9) 


and in the large /i limit these couplings can be written as 


G 


BW±W 


mz(^wS2p 
' fi{M2-Ml) ■ 


( 2 . 10 ) 


The interaction is depicted in the bottom-left graph of figure 1 . The coupling is proportional 
to and originates from the bino-wino mixing, Zi 2 . 

We illustrate the wino-bino interaction with the Z boson in the bottom-right graph of 
figure 1. Unlike the other interactions, this coupling requires both the wino-higgsino and 
the bino-higgsino interaction. Formally, the interaction is defined by the x^-Xi-Z coupling 

G^Z = ~^ijZ = - {ZisZjs — Zi4Zj4) , ( 2 . 11 ) 

S 2 W 

and reduces in the large /r limit to 


rR = ~ ^ 1 ^ 2/3 

W°BZ W°BZ ~ 2/^2 


( 2 . 12 ) 


These couplings are proportional to {mz/fJ-)‘^ and become small quicker than the other 
couplings as \fi\ increases. There exists another contribution to the bino-wino-Z interac¬ 
tion from higher order terms. The higgsino-Higgs loop diagram generates a dimension-5 
operator, 

51U“c7'^^F“ . (2.13) 
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A naive dimensional analysis suggests that the coefficient of this operator is proportional 
to al{^'K^). Although this contribution is suppressed by the loop factor, it can easily 
dominate the tree-level interaction, eq. (2.12), for very large fi. We leave a detailed study 
of effects of this operator for future work. 


3 Two-body decays 

3.1 Wino NLSP case (jMsI > |Mi|) 

Throughout this paper, we assume that gluinos are heavier than wino and bino. In this 
section we consider the cases where winos are heavier than bino. In SUSY models with 
heavy scalars and higgsinos, the mass difference between the charged and neutral winos 
is small and the decays among the wino multiplet (e.g. —)• -k^W^) can be neglected 

compared to the wino decays into a bino. We therefore do not consider the decays within 
the wino multiplet in this section. We will mention the effect of this decay mode in section 6. 

The decay mode of the charged wino is shown in the bottom-left diagram of hgure 1. 
The final state consists of and B. The decay rate of this process is given by [29]: 


T(W^ W^B) 




^ M 2 G) 

asl^M2 m| 

4 (M 2 - Mi)2 



(3.1) 


where = Mf/M^ and nw = 'niw/^ 2 - analytic expression for f-{x,y) is given in 
appendix A. Throughout this and the following sections, the first approximation is valid 
in the large y, limit, while the following assumes both fj, ^ mz and M 2 mw- Using the 
last expression, the lifetime of the charged wino is found to be 


2.5cm. (j^,^)"(5^TA)(T) ( 


M 2 


M2+M1 


(3.2) 


In figure 2(a) we show the decay rate of as a function of |/i|. The horizontal green 
line denotes the decay rate of B^ meson as a reference, where the displacement of the 
decay from the primary vertex starts to be visible. The horizontal blue line corresponds to 
the decay rate where the c/T = 1 cm, where the charged winos start to reach the trackers 
leaving the distinctive kink-like signature.^ As can be seen, the charged wino has a collider 
relevant lifetime for |/r| > 0(10^) TeV, if the two body —)■ W^B is kinematically 

allowed. 

There are two possible decay modes for the neutral wino: ZB and —)> hB. 

The former is depicted in the bottom-right diagram in hgure 1 and can be mediated either 
by the coupling eq. (2.11) or by the dimension-5 operator, eq. (2.13). However, these 

®The trackers inner radii are for example: 44 mm for CMS [33], 50.5 mm for ATLAS [34], and 16 mm 
for ILD at the planned International Linear Collider [35]. 
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Figure 2. The two-body decay rate of winos as a function of /i, with other parameters fixed at: 
M 2 = 800 GeV, Ml = 300 GeV, tan/3 = 2. 


contributions are suppressed by the extra factor and the loop factor in the matrix 

element, respectively, compared to the —)• hB decay. Therefore, if the mass difference 

is large enough to allow the two body —)• hB decay, the neutral wino predominantly 

decays into h and B. The decay rate is given by: 


T{W° hB) = 


M2'>mh 


^ Hh) 

aM2 2 

A '^ 2/3 2 Jh\hB^ hh) 


(3.3) 


where the analytic expression of fh{x,y) is given in appendix A. We 

approximate the exact formula for fi /§> mz and, additionally, M 2 /?> ruh in the second and 
third line, respectively. It is worth noting that the last expression is identical to the one 
in eq. (3.1). This can be also seen in figure 2(b), where we show the decay rate of as 
a function of |/i|. The approximate lifetime formula for the is therefore the same as in 
eq. (3.2) with cr^o ~ At colliders, the long-lived neutral wino leaves the displaced 

dijet/jets signatures [36, 37]. 

Finally, we would like to comment on a special case when mz < \M2 \ — \Mi\ < m^. In 
this case the —)• hB decay is kinematically forbidden and the IT*^ decays predominantly 

into Z and B, through the bottom-right diagram in figure 1 and the dimension-5 operator, 
eq. (2.13). Since the tree-level contribution is suppressed by the extra l//i factor, the con¬ 
tribution from the dimension-5 operator dominates in the large |/r| region. This parameter 
region is interesting, because the neutral wino can be long lived with much smaller \y\. We 
leave this question for future work. 
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Figure 3. The two-body decay rate of bino as a function of /x, with other parameters fixed at: 
Ml = 800 GeV, Ma = 300 GeV, tan/3 = 2. 


3.2 Bino NLSP case (|Mi| > [MsD 

If bino is heavier than winos, the bino can decay either to i? —)• W^W^, or B —?• hW^, 
or B ^ ZW°. As we discussed in the previous subsection, the interaction B-W-Z is 
suppressed compared to B-W-W and B-W-h, and the bino decay is dominated by the W 
and h decay modes. The partial decay rates for these modes are given by: 

a 




rs«/ 


052 / 3-^1 


m| 


/-(Mw: Mw) 


/x2 (Ml - M2)2 


and 


Mi^mw a 


a 


(Ml — M2 )s|^ 


?7r| 

/x2 


1 + 


M 2 

Ml 


(3.4) 


r(B ^ hW^) = /X,) 


/x>mz 


Oi 


aMi 2 ml f r i 

A '^ 2/3 2 

^ /X 

2 


Tij-N 2 mjz M2\ 

-(Mi-M2)s2^^(^l + ^j 


(3.5) 


where /x^ = /xw = m^/Mf and Hh = m^/Mf. The approximate formulae are 

for successively applied Mi /§> mw and Mi /§> nih, respectively. If Mi /§> ruh, we find that 
the branching ratios for different decay modes are approximately equal: 


BR{B IT+IT“) = BR{B IT“1T+) ~ BR{B hW^) ~ 1/3. 


(3.6) 


Figure 3 shows the total and partial decay rates for the bino as a function of |/x|. One 
can see that the bino can be long-lived on collider time-scales if |/x| > 0(10®) TeV. The 
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f'L 



L/R 


fL/Fi 


Figure 4. Diagram for —>■ f fB three-body decay. 


bino lifetime is approximately given as 

cr- ~ 1cm . ( ^ (±Y ( V 

^ Vl06 TeV/ \Mi + M 2 j 


(3.7) 


4 Three-body decays 

In this section we consider the cases where the mass difference between bino and wino is 
small and two-body decay modes considered in the previous section are kinematically not 
allowed. 

4.1 Wino NLSP case (jMs] > |Mi|) 

We start with a case of IM 2 I > |Mi| and define the mass difference AM = IM 2 I — |Mi| > 0. 
If AM < niw the charged wino two-body decay, —)• W^B, is forbidden. In this case, 
the charged wino decays into a pair of fermions and a bino, via an off-shell W as shown in 
figure 4. The decay rate of this process is given by: 

Tf{W+^ fJ'B) = 

li'>raz 


AM<M 2 


where jjlw = 'mw /^2 an analytic 
eqs. (A.5)-(A.9). The last approximation assumes AM <C M 2 together with /r ^ rnz- 
One may naively expect that the decay rate should be proportional to {AM)^/m'^ due to 
the s-channel structure. However the wino-bino mixing, Z 12 , is inversely proportional to 
AM, see eq. (2.5), and the decay rate is finally proportional to (AM)^. 

We show the decay rate, T{W^ —^ ff'B) = ffB), in figure 5(a) 

and 5(b) for AM = 50 and 20 GeV, respectively. It is clear that the charged wino can be 


a‘^M2 1 ^2 o ^ N 

IOtt 

a^M2 8% m| m^H_(/Xg,/xw) 

IOtt - Mi)2 

W{AMfs^ 

157r IJ? ^ ^ 

expression of Vt,±{x,y) is given in appendix A, 
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Figure 5. The three-body decay rate of the charged wino as a function of |/j,| for two values of 
AM = 50 GeV and 20 GeV, left and right panel, respectively. The other parameters are fixed at: 
M 2 = Mi + AM, Ml = 500 GeV, tan/3 = 2. 


long lived for |/i| > 0(10^ ‘^) TeV depending on the mass splitting AM. The approximate 
formula for the charged wino lifetime is given by: 


cr, 


w± 


~ 1 cm' 


10^ TeV 


30 GeV\3/ 1 n2 
AM ) 


(4.2) 


We now turn to the three-body decay of the neutral wino for AM < mz- In the 
previous section we have noted that the W-B-Z coupling is suppressed by an additional 
mz/fJ- factor compared to the W-B-h coupling, and therefore the decay mode to the Higgs 
dominates in the neutral wino two-body decay. In three-body decay, however, the Higgs 
exchange diagram also receives a suppression, which is proportional to the mass of the final 
state fermions as shown in the upper diagrams in figure 6. 

There are several competing contributions depending on the value of |/r|. Firstly, there 
is a Z-exchange diagram shown in the lower-left panel of figure 6. The matrix element 
of this diagram has the LL and RR chirality structure for the final state fermions and is 
proportional to l//r^. Additionally, if sfermions masses are of the same order as |//|, the 
sfermion exchange diagram shown in the lower-right panel of figure 6 can provide a sizeable 
contribution. The matrix element is inversely proportional to the squared mass of the left- 
handed sfermions, m?- . The contribution can interfere with the Z-exchange diagram since 

JL 

the matrix element has the LL structure. These contributions become significant if |/x| or 
m‘\ are relatively small. 

JL 

Let us assume that sfermions and higgsinos have the same mass scale, Mx, which is 
much larger than mz- Neglecting the contribution induced from the dimension-5 operator, 
the total decay rate can be written as 


r(w° ^ fjB) = r(2) (wo ^ ffB) + (it° ^ ffB) , (4.3) 
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Figure 6. Diagrams contributing to —)■ ffB three-body decay. 


where is the Higgs exchange contribution which scales as whilst 

r(‘) = y, (f f + f f + ry) is the contributions from the Z-exchange, sfermion exchange 
and their interference, scaling as 1/M^.^ 

The analytic expression for the Higgs exchange contribution, T^, is given by: 




a^M2 1 


327r 

j 


’W 


HILJ 
mw ^ 


^W°Bh ' 


fj.,mA'>mz a^M2 ^213 , 

/r/j) 


Stt 


^ 2 W 


am<M 2 4a^ (AM)' 


’ 2/3 


ISvr 


ml 


_ 2 

„2 2 
° 2 W F 


(4.4) 


where = c^j and = —Sajc^^ both of which are reduced to —1 in 

the decoupling limit of the SUSY Higgs bosons {jua 3> m^). An analytic form of yih{x,y) 
is given in appendix A, eqs. (A.lO)-(A.ll). In the last step we approximate the expression 
assuming AM <C M 2 . 

^ The .Z-exchange contribution induced from the dimension-5 operator has scaling 1/Mx in the matrix 
element and interferes with the tree-level Z and sfermion contributions. We neglect this contribution in this 
study because the operator should be suppressed by the loop factor O^ajA^n), which is seemingly smaller 
than the mf/mw factor for f — b,T in the tree-level Higgs exchange contribution. The full calculation 
including the higher order .Z-exchange contribution is however important to precisely determine the lifetime 
of the neutral wino, which is beyond the scope of this paper. 


10 







In the limit of heavy sfermions and higgsinos, fj, f ^ and are given by 




ty{W^ ^ fjB) 




ryw^^ffB) 


32tt 2si, 


hw ^ 




am<M2 f C2/3 (AM)^ 


l ,, mj:^:$>mz a ‘^ M 2^2 ^''^2 


Ml 


327r 


-K 


7 r 




Jl 

am<M2 


-/tf ■ 

127r ^ 


m 


7 l 






327r 
„2 


S2V1/ //2 

JL 


am<M 2 ^ ^ C 2/3 (AM)® 

- -2—7' 

loTT S2 m/ ml n"' 


(4.5) 


(4.6) 


(4.7) 


with liz = m^/M|, kf = |e{^||e{ 2 l and Vf = 4(13 — e/s^), where Ig and e/ are the 
isospin and the electric charge of the sfermion /. The represents the coupling between 
the left-handed sfermion and the neutralino y?, which can be written as 




jf — pfo ‘2 

l -'3 

cwsw 



(4.8) 


The y is given by (2/3 — 4 eys^)^ -|- (2/3)^, namely for each fermion we have 

u : (1 - 2'Sw)^ + 1 

d : (—1 -|- + 1 

u : 2 

e : (—1 -|- 4s^)^ -|- 1. (4-9) 


Analytic forms of 12+, Qf and ^vf are given in appendix A. 

The contributions > ffB) and r7)(IT0 —^ ffB) as well as the total decay 

rate r(IT° —)■ ffB) are shown in figure 7(a) and 7(b) for AM = 50 and 20 GeV, re¬ 
spectively, as functions of Mx, assuming ^ = mj^ = Mx- One can see that the neutral 
wino can become long lived with Mx ^ 20 TeV for the AM = 50 GeV case and even for 
Mx ^ 5 TeV for the AM = 20 GeV case. The contribution from the Z and sfermion ex¬ 
change diagrams, rO)^ is negligible for very large /r. As Mx is decreasing, this contribution 
starts to be significant around Mx ~ 20 TeV and becomes the leading contribution from 
Mx<8 TeV. 
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(a) 


(b) 


Figure 7. The three-body decay rate of the neutral wino as a function of Mx{= |/r| = for 
two values of AM = 50 GeV and 20 GeV, left and right panel, respectively. The other parameters 
are taken at M 2 = Mi + AM, Mi = 500 GeV, tan/3 = 2. 



u 

d 

e 

V 

Kj 

0.16 

0.16 

1.44 

1.44 

kfVf 

0.56 

0.68 

-1.33 

-2.40 


1.16 

1.49 

1.01 

2 


Table 1. Approximate values of coefficients kj, kfVf and appearing in eqs. (4.5)-(4.7). 

In the limit of mj^,mA 3> mz and AM ^ M 2 , the lifetime of the neutral wino can 
be written approximately as 


cr. 


wo 


>20/ sin(2/3) TeV 


1 cm • 


^ f ( 

MO GeV' 

)V> 

(4.10) 

VIO 2 

TeV/ ' 

i AM / 

1 cm • 

/ 


MO GeV' 

]'( ^ V 

(4.11) 

VIO 2 

TeV; ' 

i AM i 



with 






/i \2/1.56TeV\2 ^ ^ /1.53TeV\2- 

m? J V m-f / “ *^2/3 Mi Mi ( ^ ) • 

jLi JL^ JL^ 

(4.12) 


For ||u| mj^ > 20/ sin(2/3) TeV, the total decay rate is dominated^ by F^^^ and in this 
regime, the lifetime of the neutral wino can be approximated by eq. (4.11). Note that 
Kj can be 0(1), table 1, and in case when there is an additional large hierarchy between /x 
and the third term of eq. (4.12) might again become relevant. 


At |/x| = rrif — 20TeV and for tan/3 = 2 we have = 0.195. 
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Figure 8 . The three-body decay rate of the bino as a function of |/r| for two values of AM = 
50 GeV and 20 GeV, left and right panel, respectively. The other parameters are taken at Mi = 
M 2 -f AM, M 2 = 500 GeV, tan/3 = 2. 


4.2 Bino NLSP case (|Mi| > IM 2 I) 

In this subsection we assume bino is heavier than winos and define AM = |Mi| — IM 2 I > 0. 
If AM < mw, only three-body bino decays are allowed. As discussed in the previous sub¬ 
section, the Higgs exchange diagram is suppressed by the mass of the final state fermions, 
and the IV-exchange diagram for B —)• ffW^ dominates bino decay. The partial decay 
rate can be written as 

TfiB^ffW^) = 


AM<Mi 


levT 

Ml m?2; '^w^- ) 

levr IJ? {Ml - M2Y 

4a^ (AMf 4/3 
IStt ’ 


where /j,]y = m^/M^ and we approximate the expression assuming AM <C Mi on top of 
// 3> mz in the final step. An analytic formula for {fi^, fiw) is provided in appendix A, 
eqs. (A.5)-(A.9). 

In figures 8(a) and 8(b) we show the bino decay rate r(.B — )• f f'W4 = J2f ^ f{B —t 
ff'W4 cts a function of |/r| for AM = 50 and 20 GeV, respectively. As can be seen, the 
bino becomes long lived for \fj,\ > 0(10^“^) TeV depending on AM. The approximate 
formula for the bino lifetime is given by 


CT^ ~ 0.5 cm ' 




104 TeV 


2/30 GeV\3 


AM 


(-)' 

^82/3^ 


(4.14) 


Note that if the decay is mainly mediated by the Higgs boson exchange, the final state will 
be dominated by (possibly displaced) 6 or r jets. 
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5 Large mass splitting between gauginos and higgsinos 


In the previous sections, we found that wino and bino can be long-lived for \fi\ > 0(10^“®) 
TeV when two-body decays are allowed, and |/x| > 0(10^“'^) TeV for three-body decays. On 
the other hand, the wino and bino masses should be less than 0(1) TeV to be produced at 
colliders. The long-lived wino and bino in this scenario can be observed only if an enormous 
mass splitting between gauginos and higgsinos is realised. 

The higgsino mass, /r, is the only dimensionfull parameter in the MSSM Lagrangian. 
Since it is a supersymmetric mass term, its origin can be different from the soft masses of 
gauginos and scalars. If one does not ask about the origin of the /i term and does accept the 
hne tuning in the electroweak symmetry breaking, phenomenologically there is no problem 
in taking at any value, as long as phenomenologically correct EWSB is achieved. 

If ^ is very large, the low energy effective lagrangian is obtained by integrating out the 
heavy higgsino fields. In doing so, the electroweak gauginos receive the threshold correction 
from the heavy higgsinos as [38, 39] 


AMf = —Il 


Att 5 


AM2^ = — L 

dvr 


(5.1) 


with 


L = 


l/iP — ' ’ 


(5.2) 


where uia is the CP-odd Higgs mass. If ruA and jjL are of the same order and j^j > 0(10^) TeV, 
these threshold corrections exceed 0(10) TeV and the Wino and Bino become out of the 
collider reach. In order for these corrections to be small and realise the large mass splitting 
between heavy higgsinos and light gauginos, one needs ruA |a*|- Since ruA is related to 
IX by 

-I-, (5.3) 


some unnatural cancellation may be required to accommodate observable long-lived winos 
and binos for |/i| > 0(10^) TeV. 

Another possibility is to attribute the origin of /r to the SUSY breaking. For example, 
in gravity mediation, /x can be generated as [40, 41] 


^ = cm3/2, (5.4) 

where c is a dimensionless coefficient. In this scenario the gaugino masses receive the 
anomaly mediated contribution [38, 39], 

= ^"i3/2 , = -^3m3/2 , (5.5) 

as well as the higgsino threshold correction of eq. (5.1). The size of these correction exceeds 
0(10) TeV for |/x| > O(IO^) TeV. In this scenario, some unnatural cancellation among these 
contributions may be required for the long-lived bino and winos within the collider reach 
for \n\ > O(103) TeV. 
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From these considerations one may find the case of the wino NLSP with IM 2 I — 

I Ml I < mz particularly interesting because the neutral wino can become long-lived for 
1^1 > 0(10^“^) TeV.® With |/i| ~ 0(10^“^) TeV, the contributions from the higgsino 
threshold corrections and the anomaly mediation are of the order of 0(10^“^) GeV and a 
large variety of the wino and bino spectra, including compressed spectrum for the three- 
body decays, can be naturally achieved by the interplay between these contributions [7, 43- 
45]. In order to precisely predict the lifetime of the neutral wino in this regime, the calcu¬ 
lation of the Z-exchange diagram induced by the dimension-5 operator may be necessary. 
We leave this task for future work. 


6 Collider signatures 


In supersymmetry with heavy scalars and higgsinos, the bino cross section is largely sup¬ 
pressed both at hadron and e'’"e“ colliders. The binos may nevertheless be produced from 
decays of gluinos at hadron colliders. As shown in sections 3 and 4, the lifetime of bino 
can become larger than 0(1) cm for |^| > 0(10®) TeV if the two-body decay is allowed 
and \fi\ > 0(10^) TeV otherwise. 

If squarks are heavy the gluinos may also become long-lived. The gluino lifetime is 
given by [8]: 


1 TeV\5/ rriq 




, 103 TeV 




( 6 . 1 ) 


Therefore, if the squark masses and |/r| are of similar order, the long-lived binos are pro¬ 
duced from long-lived gluinos. In order to have the gluino decay well inside the detector, 
|/i| and the squark mass cannot be too large. A number of searches have been performed 
by the LHC experiments looking for metastable gluinos that hadronizes into a colorless R- 
hadron. These include cases when the i?-hadron escapes a detector before decaying [46, 47], 
the i?-hadron is stopped in a detector [26, 48, 49], or it decays inside the detector [50]. 

Bino can become long-lived already with |/r| rsj O(IO^) TeV. In this case the bino pre¬ 
dominantly decays into an off-shell W and a charged wino. The charged wino subsequently 
decays into an off-shell W and the neutral wino, with potentially long lifetime [51, 52]. 
The charged wino lifetime in this region is given by 


^ /160MeV\3,^ 

cr^± ~ 5 cm • ( —^- 1(1- 


Am 


w 


mz 




-1/2 


( 6 . 2 ) 


where the mass splitting Am^ = — m^q can be written as 




1 


„2 ^2 




(6.3) 


at tree-level and receives a radiative correction, AmT® ~ 160 MeV [53]. The long-lived 
charged winos are searched for by looking for the disappearing track signature [54, 55] or 

® This mass spectrum of bino and wino is also motivated by a recent scan of the 10-dimensional parameter 
space of phenomenological MSSM [42]. 
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measurements of ionisation energy loss in a pixel detector [49] . If the long-lived binos carry 
the charged winos in the middle of the tracking system (50-100 mm), the signal can be 
seen as an appearing-and-disappearing track signature in the detector. 

The long-lived wino can be produced either directly at hadron and e'’"e“ colliders or 
indirectly from prompt/non-prompt decays of gluinos at hadron colliders. The charged 
wino predominantly decays into an on- or off-shell W and the LSP bino. If the decay 
products of IT=*=(*) are reconstructed, the signal may be detected as a kinked tracks in the 
leptonic channel or a displaced dijet with a charged track pointing to the secondary 
vertex in the hadronic W^*'^ channel. 

The neutral wino predominantly decays into a and a bino in most cases. The signal 
should be detected as the displaced jets/dijet signature [36, 56, 57]. It is worth noting that 
the Higgs mediation can be pinned down by confirming that the displaced jets originate 
from 6-quarks and r-leptons [58-61], as opposed to the Z mediation case, where light jets 
would dominate the final state. In the region where the neutral wino is long-lived with 
j/rj ~ 0(10^“^) TeV, the two-body decay is kinematically forbidden and the decays 
into a bino and an off-shell Z, h or left-handed sfermions with a small bino-wino mass 
splitting. Because the decay products in these off-shell decays are soft, the e''"e“ may offer 
the best opportunity to detect the long-lived neutral wino in this parameter region. At 
hadron colliders, the search for the displaced Z^*') —)> ££ may also be promising [57, 62, 63]. 

7 Conclusions 

We investigated a possibility of having long-lived binos and winos in SUSY models with 
heavy scalars and higgsinos. In the parameter region of interest the SU{2)l and U(l)y 
gaugino sectors are decoupled from each other with very small mixings proportional to n//x. 
In this region, the heavier of bino and wino practically does not interact with the lighter 
one and its lifetime becomes relevant for collider experiments, cr > 0 ( 1 ) cm. 

We revisited the decays of bino and winos and found simple formulae for the decay 
rates and lifetimes, which are valid when the scalars and the higgsinos are much heav¬ 
ier than the gauginos. We have found that the long-lived bino and wino emerge when 
\^\ > 0(10®“®) TeV if the two-body decay mode is open and ]//] > 0(10^“^) TeV other¬ 
wise. One exception is the case with 0 < jM 2 l — jMij < mz, where the lifetime of the 
neutral wino becomes 0(1) cm for j/rj ~ 0(10^“^) TeV, depending on the wino-bino mass 
splitting, because the Higgs exchange diagram is suppressed by the mass of the final state 
fermions. We compared our results to SDecay finding a good agreement for the range of 
parameters permitted by SDecay. 

We briefly discussed how the large mass splitting between gauginos and higgsinos can 
be achieved. If the origin of /r is independent of the SUSY breaking, the large mass splitting 
can be realised relatively easily, although the threshold correction from the heavy higgsinos 
to the gaugino masses needs to be suppressed. On the other hand, if the ^ is linked to the 
SUSY breaking and in particular ~ ^ 3 / 2 , the contributions from the higgsino threshold 
correction and the anomaly mediation become significant. However the large splitting is 
still possible if one arranges the cancellation between these contributions. 
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We also discussed a possible collider signature for the long-lived bino and wino in 
this scenario. The production of bino is only possible from the decay of gluinos, although 
gluinos tend also to be long-lived for mq > O(IO^) TeV. For the long-lived bino NLSP 
case, the charged wino may also be long-lived, because of the very small mass splitting 
within the wino multiplet. In S —)■ decay, the bino may carry the long-lived 

charged winos into the middle of the tracking system and the signal could be seen as an 
appearing-and-disappearing track signature. 

For the wino NLSP case, the production of winos is possible either directly or indirectly 
from the gluino decay. The long-lived charged wino decaying to and B can be 

detected as a kinked-track signature from the leptonic decay of W^*\ The hadronic 
mode may also be seen as events with a displaced dijet and a track pointing to the secondary 
vertex. Detecting the long-lived neutral wino with \ij,\ ~ 0(10^“^) TeV may be challenging 
because the decay products will be soft due to the small wino-bino mass splitting. In 
this case, e'*'e“ collider may be ideal to detect the long-lived neutral wino. Otherwise a 
displaced off-shell Z* —)• ii decay may be promising even at hadron colliders. 
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A Auxiliary functions 

In this appendix we summarize the analytic expressions for functions used in the calculation 


of decay widths. 

f±{x,y) = VX{x,y)ri±{x,y) , (A.l) 

A(a:, y) = 1 + + y^ - 2x - 2y - 2xy , (A.2) 

V±ix,y) = (1 +X - y) + — —— - — —±6^/x:, (A.3) 

y 

fh{x,y) = {l + x-y + 2y/x)^/M^x^, (A.4) 

^±{x,y) = 2F{x,y)±G{x,y), (A.5) 

X — 1 

F{x, y) = [w(x, y) -f y(5 -f 5x - ly)] 

-|{(1 + x-y)h\x + [a;(x, y) + 2x]£(x, y)} , (A.6) 


G{x, y) = ^/x[A{x - 1) + {I + x - 2y)h\x + {uj{x,y) - y{l + x - y)}C{x,y)\ , (A.7) 
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a;(x, y) = 1 - 2x - 2y + (y - x)^, 


(A.8) 

2 


^ 1 - X - y 

, (A.9) 

y/-U}{x,y) 

L \^-u{x,y)J 

^y/-a;(x,y)^J 


^h{x,y) = H{x,y) + G{x,y), 


(A.IO) 


H{x, y) = ^(1 — x){6y — 5 — 5x) + ^(1 — 4y — 4xy + 3y^ + x^) Inx 

+ ^(-5x^y - 3y^ + 7y^ + 1 - - a: + - 5y + 7xy^ - 2xy)£(x, y ), 

(A.ll) 


Qf{x) 


^(1 — 8x + 8x^ — x^ — 12x^ Inx) 

+ Tl + x(9 — 9x — x^) + 6x(l + x) Inxl , 
6 


Qv{x) = —f — 5 — 9x + 9x^ + 5x^ + (3 — 9x — 18x^) Inx) 
36 V / 


(x — l)yjAx — 1 


arctan 


\/4x — 1 


+ arctan 


1 - 2x 
\/4x — 1 


(A.12) 


(A.13) 


^vf{x,y) = (1 + \/x) 

[l - 3x - 3x^ + x^ - 3(1 + x)y^ + 2y^] 

+ [5 + 5x2 _ ^ 2 (y _ i)y + 2x(l + 6y)] ^ 

+ [l + x2 + x(y - 2) + y - 2y2] [x^ + (y - l)^ - 2x(l + y)] 

^_ M{x,y) _ 1 

6 y^-x2 - (y - 1)2 + 2x(l + y)-l ’ 


(A.14) 


M.{x,y) = arctan 


-1 + X + y 


— arctan 


Ly^-x^- (y-l)2 + 2x(l + y)J 
1 - X + y 


V^-x2 - (y - 1)2 + 2x(l + y)J 


(A.15) 


B Comparison with SDecay 

In this section we compare approximated results obtained in this paper with exact calcu¬ 
lation performed using the computer code SDecay [27]. Figure 9 shows the partial decay 
widths using both approaches for processes: —)■ W^B (cf. section 3.1), W~^ —)■ e^vB 
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(cf. section 4.1), and —?■ bbB (cf. section 4.1). For the last process we remove the squark 

mediated contribution by taking = 10|//| and approximate F = + Ff), where 

Nc{= 3) is the colour factor. We use eqs. (4.4) and (4.5) for F^ and F^, respectively. In 
addition, we show the ratio of decay widths, 

■papprox 

D _ _ _ 

pSDecay ’ 

in the lower panels of each plot. In order to ensure that the spectra and couplings are the 
same for the purpose of comparison, we calculate the spectrum fed to SDecay externally 
to avoid corrections due to running parameters. For sufficiently large values of /r we find 
a very good agreement between both predictions, within 0(10%). A level of 0(40%) 
agreement can be also achieved for lower values of /i. An important thing to note is that 
the calculation by SDecay fails for /i > 10® TeV, due to numerical problems. Another point 
to note is that because of the large separation between gaugino and stops scales one would 
typically run into problems with standard spectrum calculators. For these reasons using 
approximate formulae would be preferred over automated codes in the large /r regime. 
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